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Abstract 



Undulatory locomotion is common to nematodes as well as to limb- 
less vertebrates, but its control is not understood in spite of the iden- 
tification of hundred of genes involved in Caenorhabditis elegans lo- 
comotion. To reveal the mechanisms of nematode undulatory loco- 
motion, we quantitatively analyzed the movement of C. elegans with 
genetic perturbations to neurons, muscles, and skeleton (cuticle). We 
also compared locomotion of different Caenorhabditis species. We con- 
structed a theoretical model that combines mechanics and biophysics, 
and that is constrained by the observations of propulsion and mus- 
cular velocities, as well as wavelength and amplitude of undulations. 
Wc find that normalized wavelength is a conserved quantity among 
wild-type C. elegans individuals, across mutants, and across different 
species. The velocity of forward propulsion scales linearly with the 
velocity of the muscular wave and the corresponding slope is also a 
conserved quantity and almost optimal; the exceptions are in some 
mutants affecting cuticle structure. In theoretical terms, the optimal- 
ity of the slope is equivalent to the exact balance between muscular 
and visco-elastic body reaction bending moments. We find that the 
amplitude and frequency of undulations are inversely correlated and 
provide a theoretical explanation for this fact. These experimental 
results are valid both for young adults and for all larval stages of wild 
type C. elegans. In particular, during development, the amplitude 
scales linearly with the wavelength, consistent with our theory. We 
also investigated the influence of substrate firmness on motion pa- 
rameters, and found that it does not affect the above invariants. In 
general, our biomechanical model can explain the observed robustness 
of the mechanisms controlling nematode undulatory locomotion. 
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Introduction 

Undulatory locomotion enables invertebrates such as nematode worms and leeches, 
and limbless vertebrates such as snakes and fish to move through different environments 
including soil, sand, water, and tissues of plants in search of food (Gray, 1964). Nema- 
tode body undulations are controlled by the neuromuscular system, which generates a 
wave of longitudinal muscle contractions modulated by elastic properties of the cuti- 
cle and hydrostatic skeleton. Such a wave consists of alternating phases of dorsal and 
ventral muscle contractions (worms lie on their sides) that travel posteriorly along the 
body length if the worm moves forward, and anteriorly if the worm moves backward. 
As a result of the interactions between neurons, muscles, skeleton/cuticle, and the en- 
vironment, the worm crawls and its body follows approximately sinusoidal trajectory 
(Fig. 1). 

There remains a major gulf between our understanding of biomechanics (Alexander 
and Goldspink, 1977; Gray and Lissmann, 1964; Niebur and Erdos, 1991; Wu, 1971; 
Cheng et al, 1998; Hirose, 1993; Ijspeert, 2001) and molecular genetics (Bargmann, 
1998; Francis et al, 2003) of undulatory movement. Studies of leech and lamprey 
(Friesen and Gang, 2001; Skinner and Mulloney, 1998; Lockery and Sejnowski, 1992; 
Marder and Galabrese, 1996; Williams, 1998; Gortez et al, 2004) have led to a system- 
level, mostly neural, understanding of such movement but have little connection to how 
such movement is specified by the genome. Molecular genetic studies have identified 
hundreds of genes involved in locomotory behavior in C. elegans and while functional 
connections among many of these genes have been elucidated, they have not yet ex- 
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plained how the system works. Sensory behavior has begun to yield to a molecular 
approach (de Bono and Maricq, 2005), but remains a long way from motor output. 
We seek to understand how coordinated motor output is specified by a genome, and 
specifically, how genes infiuence the parameters that directly control the locomotion. 
We therefore started a two-pronged approach to this general problem in C. elegans. 
First, we started with single gene perturbations and measured behavior quantitatively. 
Second, we started building relatively simple models of worm movement, matching 
main parameters to experimental observables. This interactive approach allowed us to 
construct a biomechanical model that fits experimental data on the characteristics of 
sinusoidal locomotion. 

The experimental approach consisted of genetic perturbations to parameters rele- 
vant for locomotory control in C. elegans that included neurons, muscles, and cuticle, 
and quantitative data analysis of the resulting motion. We also compared locomotion 
of several wild-type Caenorhabditis nematode species to provide additional clues about 
generic characteristics. Our theoretical model combines mechanical and biophysical 
aspects of undulatory locomotion. This model allows us to explain why some charac- 
teristics change or stay constant as a genetic perturbation is applied. In this respect, 
our model has a potential to provide richer information than standard, neural-level ap- 
proaches in leech and lamprey (Priesen and Gang, 2001; Skinner and MuUoney, 1998; 
Lockery and Sejnowski, 1992; Marder and Calabrese, 1996; Williams, 1998). 

We find that although several quantities characterizing the movement vary from 
mutant to mutant, their inter-relationships are preserved in most cases, suggesting that 
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locomotory control mechanisms are robust and somehow evolutionarily optimized. In 
particular, we find that our theory is consistent with the experimental findings that 
(i) the velocity of forward locomotion scales linearly with the velocity of the muscular 
contractions wave such that the former is close to optimal for all species and almost all 
mutants we examined, (ii) the wavelength of undulations, normalized with respect to 
body length, is highly conserved across different related species and different mutants, 
and (iii) the amplitude of the undulations weakly decreases with the frequency of the 
wave. These results are preserved during different developmental stages of C. elegans, 
as well as on different substrate firmness. 

The system responsible for locomotion of Caenorhabditis worms comprises four main 
elements: nerve ring (head) neurons, motor neurons, complexes of muscles with skele- 
ton/cuticle, and different mechano-sensory feedback loops that influence activity of 
motor neurons (Fig. 2). The precise interactions between these elements in producing 
oscillatory locomotory output remain unknown. In this paper, we simply assume that 
this system is capable of producing an oscillatory wave of muscle contractions that prop- 
agates along the worm's body. The main results and conclusions below are independent 
of the particular oscillatory mechanism. Construction of the model is described in the 
Theoretical Model section. 
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Experimental procedures 

Mutant and species selection. 

We examined the locomotion of several C. elegans mutant classes with defects in neural, 
muscular, and cuticle functions. For mutations affecting neurons we studied: MT2426 
goa-l(nll34), PS1762 goa-l(syl92), and PS4498 egl-30(tg26), which alter G-protein al- 
pha subunit genes responsible for regulating synaptic transmission (Mendel et al, 1995; 
Segalat et al 1995; Moghal et al, 2003), cat-2(elll2), which encodes a tyrosine hy- 
droxylase, an enzyme required for biosynthesis of dopamine (Lints and Emmons, 1999) 
and cat-4(ell41), which encodes GTP cyclohydrolase I, an enzyme required in the pro- 
cess of synthesis of dopamine and serotonin (Kapatos et al, 1999). cat-2 and cat-4 
affect activities of dopaminergic neurons (Sulston et al, 1975). The specific neuronal 
mutants chosen by us differ from a vast majority of other neuronal mutants because 
of their known hyperactive locomotion. We want to investigate this interesting fea- 
ture more quantitatively. For mutations affecting muscles we studied several strains 
carrying mutations in the head region of body wall muscle myosin unc-54 RW130 
unc-54(stl30), RW132 unc-54(stl32), RW134 unc-54(stl34), RW135 unc-54(stl35), 
RW5008 unc-54(s95), and BC347 unc-54(s74)', these mutations are hypothesized to 
alter the contraction-relaxation cycle of the myosin-actin crossbridge formation by in- 
creasing its duration (Moerman and Fire, 1997). For mutations affecting cuticle, we 
studied two loss-of-function alleles of sqt-1: BElOl sqt-l(sclOl) and BE103 sqt-l(scl03) 
mutants, sqt-1 encodes a cuticle collagen (Kramer et al, 1988), a protein responsible 
for elastic and structural properties of the cuticle. We also studied a mutant BE109 
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of an unknown gene that also affects cuticle by removing the struts that separate lay- 
ers of the cuticle (J. Kramer, pers. comm.), and two mutants with increased body 
length: CB185 lon-l(el85) and double mutant PS697 lon-l(el85); lon-2(e678). The 
lon-1 gene encodes a protein belonging to the PR-protein superfamily that regulates 
polyploidization and body length (Brenner, 1974; Maduzia et al 2002). Moreover, we 
examined several wild-type nematode species closely related to N2 C. elegans (Brenner, 
1974): SB339 C. japonica (Kiontke et al, 2002), AF16 C. briggsae, SB146 C. remanei 
(A. Fodor, pers. comm.), CB5161 Caenorhabditis sp. (A. Fodor, pers. comm.), and 
PSIOIO Caenorhabditis sp (R. Giblin-Davis, pers. comm.). 

Description of experimental setup. 

We video recorded and digitized the motion of young adult hermaphrodite (or female 
from malc-fcmalc strains) Caenorhabditis worms: wild-type C. elegans and their mu- 
tants, related species, and additionally C. elegans larvae. Young adults were 15-20 hr 
post mid-L4 developmental stage, and larvae were in all stages from LI to young adults. 
Adult C. elegans and related species are tiny animals with length about 1 mm and width 
about 0.1 mm, so our experimental set-up involved a microscope connected to a video 
and a specialized open source software. Initially, an agar plate was covered with thin 
film of E. coli OP50 bacteria mixed with LB media. After that, the agar plate was 
covered with water (about 0.1 mm of hight) and left for 1 hour so that the mixture of 
bacteria with LB media is absorbed in the agar and the surface dry. Thus, during the 
recording sessions worms moved not through water but through a layer of bacteria with 
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some remaining LB media. That type of movement can be classified as crawHng, since 
the worms touch the substrate (agar), and are only partially immersed in the bacteria 
layer. 

In experiments with different substrate firmness, the concentration of agar in the 
substrate solution was changed. Higher agar concentration corresponds to a more stiff 
substrate on which worms move. 

The video recording and data extraction was done using a device specially designed 
for studying Caenorhabditis locomotion (Cronin et al, 2005). We collected 5 minutes 
of video per worm, extracting digital locomotion data from the middle 4 minutes. Such 
4 minute windows average over possible sensory influences that can vary among worms 
and thus statistically minimize the variability of external conditions. From these data we 
derived values of the velocity of forward and backward motion, frequency, wavelength, 
and amplitude of undulations. Our experimental setup allows us to measure both 
instantaneous and average values of crawling parameters. 

Theoretical model 

The theoretical part of this study involved the construction of a comprehensive 
mathematical model, from which we could compute the quantities directly measured in 
the experimental part. 
Mechanics. 

The mechanical aspect of the undulatory locomotion is modeled by assuming that 
worm's body can be treated as an active bending beam (Wu, 1971; Cheng et al, 1998). 
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Crucial for the motion generation is the existence of neuromuscular wave that propa- 
gates along worm's body. In C. elegans such a wave can be generated by long-range in- 
teractions between body segments, as it was shown theoretically for a system of coupled 
oscillators (Kopell and Ermentrout, 1988). Prom a biological perspective, long-range 
interactions could be mediated by stretch receptors, which are hypothesized (in White 
et al, 1986; Chalfie and White, 1988) to be located on endings of the extended dendritic 
processes of the excitatory motor neurons (Fig. 2). These extensions can be as long as 
25 % of the worm's length. For the forward motion B neurons, these processes extend 
posteriorly, i.e., the coupling via stretch receptors comes from the posterior parts of 
the body. For the backward motion A neurons, the extended processes are directed in 
the opposite direction, and consequently the stretch receptor coupling comes from the 
anterior parts of the body. This distinction in dendrite directionality is correlated with 
observed opposite directions of neuromuscular waves in both locomotory circuits. 

Rhythmic muscle contractions, caused by the motor neurons' activity, bend the 
worm's body and this generates the propulsion of its center of mass. Newtonian equa- 
tions of motion describe the balance between muscular, elastic, frictional, and internal 
pressure forces. During crawling, inertial effects are negligible because of the small mass 
of the worms (Gray, 1964) and small maximal acceleration (for estimation see, Niebur 
and Erdos, 1991). Equations of motion are written for a slice of a worm, perpendicular 
to its main body axis, and for a small slope of undulation. We choose a system of 
coordinates such that at resting conditions when the worm's body is straight, its main 
body axis is parallel to the x coordinate (Fig. 3). The transverse force balance is 
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F+{F + 5F) cos(50) + {G + 6G) sm{-5(p) - F^Sx = 0, 



(1) 



and the longitudinal force balance is 



G+{G + SG) cos{S(f)) - (F + SF) sm{-S(f)) + FlSx = 0, 



(2) 



where F is the beam shear force, G is the longitudinal tension (it includes both tension 
in structural elements of cuticle and internal hydrostatic pressure), 50 is an infinitisemal 
change of undulations angle over the slice width, and F^ and Fl are the normal and 
longitudinal components of the frictional force per unit length of the shce. We assume 
that these components are proportional to the normal vn and tangential components 
of the animal's velocity, with proportionality coefficients cn and cl, respectively. With 
the sign convention as in Fig. 3, we have F^ — —CNVjsf and F^ — clv^. This assumption 
is analogous to an assumption made in fluid dynamics models, and it is known as 
"Resistive Force Theory" (Gray and Hancock, 1955; Lighthill, 1976). The normal and 
longitudinal components of velocity can be represented by main body axis velocity v (it 
is parallel to the main body axis), lateral velocity u = dh/dt with h being the lateral 
displacement, and the tangential angle to the slice in the form: — v cos (p — usin 0, 
vn — vsincf) + u cos 0. 

After expanding cos(50) and sin((50) for small 50 in Eqs. (1) and (2), and neglecting 
higher order terms, we obtain the following differential equations of motion: 
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dF d<t^ 
'dx ~ dx ^ 



(3) 



and 



dG , , 

^TT + ^ = -CLVL{x,t). 
OX OX 



(4) 



The beam shear force F is related to the total bending moment M by 



F = -dM/dx. (5) 

Equations (3)-(5) constitute the basis for our analysis of mechanics of worm's un- 
dulatory locomotion. They can be further simplified in the limit of small angle of 
undulations cf). In this hmit sin0 tan0 = dh/dx and cos0 f» 1 — 0^/2, and 
as a result velocities take the forms: vl ^ v 1 — \{dh/dxY — (dh/dt) (dh/dx), and 
vn ~ vdh/dx + dh/dt. Also, the body curvature dcp/dx ~ d'^h/dx^ a; 0, and there- 
fore, the term containing longitudinal tension G in Eq. (3) can be neglected. Using 
this information, we can combine Eqs. (3) and (5) to arrive at the equation describing 
spatio-temporal pattern of the total bending moment: 



d'^M{x,t) 
dx"^ 



cn 



dh dh 

) \ 

dx dt 



(6) 



One usually solves this equation for M given some sinusoidal form of the lateral dis- 
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placement h{x,t) that mimics undulatory locomotion (Cheng et al, 1998). This is 
inverse to a seemingly more "natural" approach, in which some form of the bending 
moment would be chosen and the equation solved for h{x,t). The general problem 
with the latter approach is that it is hard to guess the right form of M{x,t) which 
produces stable undulatory motion. In particular, the spatial dependence of M is cru- 
cial for stabilization, and, as it turns out, M depends nontrivially on the position x 
along the body (see below). Also, it is not clear what boundary conditions to im- 
pose on the h{x, t) function, since during real undulatory motion the animals' head 
and tail are in permanent motion. For these reasons, we adopt the inverse approach 
of Cheng et al. (Cheng et al, 1998) and solve Eq. (6) for M with boundary condi- 
tions M(0,t) = M(L,t) = dM{0,t)/dx = dM{L,t)/dx = 0, where L is the worm's 
body length. These conditions follow from the natural requirement that all forces and 
moments must disappear outside the borders of worm's body. We choose the lateral 
displacement h in the form: 

h{x, t) = Ao cos{ujt - 2nx/X) + Sh{x, t), (7) 

where Aq is the amplitude of undulations, cu and A are the angular frequency and 
wavelength characterizing the neuromuscular wave which travels to the right with the 
velocity Aa;/27r, and 5h{x,t) — A{t) -\- B{t)x is the so-called recoil correction (Cheng 
et al, 1998). This correction ensures that the total external force acting on worm's 
body is zero at all times (Lighthill, 1960; Pedley and Hill, 1999). From a mathematical 
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point of view, this correction is necessary to satisfy the boundary conditions, and the 
functions A{t) and B{t) are determined self-consistently from them. From a biological 
perspective, one can view this correction as a mechanosensory feedback contribution, 
which is necessary for coordination of the movement. 

The total bending moment M is composed of two additive contributions: M = 
Mm + Me, where the muscle generated bending moment is given by Eq. (12) (see 
below), and visco-elastic reaction of the body moment Mg is given by: 

Me = -EJd^h/dx'' - nJ^{d^h/dx^), (8) 



where E is the Young's (elastic) modulus of the hydrostatic skeleton and cuticle, J is the 
inertial moment of the body in relation to the axis perpendicular to the body midline, /j, 
is the viscous coefficient characterizing relaxation processes in the hydrostatic skeleton 
and cuticle. 

The propulsion velocity is equivalent to the main body axis velocity v. Having the 
total bending moment M, we can find the equation for the propulsion velocity from 
Eqs. (4) and (5). First, we substitute Eq. (5) for the beam shear force F in eq. 
(4), and then we use the fact that in the limit of a small slope of undulations ~ 
V 1 — ^{dh/dxy — [dh / dt){dh / dx) . Both corrections to or w are small of the order 
of 0{h^), however only the second, {dh/dt){dh/dx), is relevant for the determination 
of V in the 0{h^) order; the first correction provides a higher order contribution 0{h^). 
This allows us to derive v from a resulting equation. However, because both M and 
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h are functions of space and time, we must perform spatial and temporal averaging of 
this equation. Spatial averaging is done within the limits of and L, where L is the 
worm's length, while temporal averaging is performed over one period of oscillations. 
After all the steps we obtain: 



where the bracket (...) denotes temporal averaging, and we neglected the contribution 
coming from the longitudinal tension, {G{L) — (7(0)), because of a partial or complete 
cancelation of {G{L)) and (G(0)) in the limit of small slope of undulations. Equation 
(9) allows us to determine self-consistently v, since the right hand side of it also depends 
on V (via M). Note that both terms under the integral in Eq. (9) are of the same order, 
i.e., 0{h?), since the total bending moment M is of the order 0{h) (see Eq. (6)). This 
implies that the propulsion velocity v is Oih?) order quantity. 

Muscle dynamics. 

Muscle activity is modeled as a simple first-order kinetic equation, in which an external 
driving factor activates muscles. That factor originates from neural oscillations at the 
neuromuscular junction and subsequent calcium influx. Since the magnitude of neural 
oscillations can in general be position dependent (as is the total bending moment M; 
see below) and we do not know in advance what that dependence could be, in what 
follows, we will compute only space averaged (over the whole worm's length) quantities 
associated with muscle activation. We use the following equation: 
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dna{t) 
dt 



(10) 



where rT^{t) is the space-averaged muscle activity level at time t, a denotes either 
dorsal (d) or ventral (v) side of the worm, Da{t) is the average driving factor, and 
the parameters / and g correspond to the rates of muscle activation and deactivation, 
respectively. We choose the space-averaged driving factor Da{t) in the form: 



where the sign -|-(— ) corresponds to dorsal (ventral) side, L is the worm's length, T 
is the constant characterizing the overall amount of calcium/synaptic transmission at 
the neuromuscular junction coming from neural oscillations, and ^ is some phase factor 
reflecting temporal delays in neural and body reaction activities. 

There are two equivalent ways of thinking about what constitutes muscle activity. 
Prom a biophysical perspective one can interpret muscle activity as the level of muscle 
contraction, which is related to the fraction of crossbridges created (Huxley, 1957). In 
this model, myosin and actin molecules interact to create crossbridges between thick 
and thin filaments. The crossbridge dynamics result in both filaments sliding past one 
another, which leads to muscle contraction and force generation. In this interpretation, 
/ is the effective rate of crossbridge association and g is its dissociation rate. Alter- 
natively, from a physiological perspective, one can think about rr^{t) as the position 



D^{t)^{T/2L) [l±cos(a;t + 0], 



(11) 
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averaged muscle membrane potential. 

Muscle dynamics and mechanics can be coupled by relating the muscle activity 
to the spatial average of the muscle bending moment Mm- We assume that the average 
muscle bending moment Mm{t) is proportional to the difference in average activities 
of the dorsal and ventral muscles with a proportionality constant k measuring 
muscle stiffness, i.e., 

M;;,{t) = K[n-d{t)-n-,{t)]. (12) 

This coupling allows us to derive the amplitude of undulations Aq as a function of 
the frequency. This can be done by balancing spatially averaged bending moments 
derived from the mechanics and muscle dynamics parts. The details are presented in 
the Appendix. 

Experimental results 

Primary locomotory data. 

In Tables 1-3 we present primary data we collected, which includes average values of the 
velocity of propulsion v, frequency of undulations a;/27r, amplitude of undulations Aq, 
wavelength of the body posture A, the ratio of the wavelength to the length of the body 
X/L, and the ratio of the propulsion velocity to the velocity of the neuromuscular wave 
denoted by 7 (the efficiency coefficient; see below). Although, many of these parameters 
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change under different conditions, tfiere are some regularities in tfie data. 

Table 1 contains locomotory data for wild-type C. elegans during its development 
through different larval stages to young adult. We observe that worms' velocity in- 
creases monotonically as they mature and this is accompanied by monotonic changes 
in amplitude and wavelength, with frequency exhibiting some variability. Contrary to 
that trend, the body-length normalized wavelength X/L and the efficiency coefficient 7 
stay relatively constant; about 2/3 for A/L, and 7 in the range 0.8 — 0.9. 

In Table 2, we present data for adult C. elegans mutants and related species, with 
wild-type C. elegans as a control. Most hyperactive mutants with affected nervous 
system move with greater velocity than wild-type. The outlier from this trend cat-4 is 
an interesting case because it moves slower but with a higher frequency as compared to 
wild-type. It is important to note that our neuronal mutants are exceptions, because 
the majority of other neuronal mutants (not examined by us) are sluggish, which is 
the opposite behavior. Mutants affecting muscle structure and properties of the cuticle 
move slower than wild-type and hyperactive neuronal mutants. The locomotory data 
show high variability in velocity, frequency, and amplitude across different mutants and 
species. However, again, the normalized wavelength X/L and the efficiency coefficient 7 
do not change much, especially X/L which is around 2/3. In most cases the coefficient 
7 is in the range 0.7 — 1.0. There are few exceptions with 7 0.5 — 0.6. These include 
cat-4, sqt-1, and one strain of unc-54- The latter is probably a consequence of the fact 
that these worms move very slow, often pausing, with mean velocity comparable to its 
standard deviation. 
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Table 3 presents locomotory data for adult C. elegans on substrates with different 
firmness (solidification) properties that can be changed by agar content in the substrate 
solution. We observe that as substrate becomes more firm, velocity and frequency 
decrease monotonically, suggesting that it is harder for worms to move forward. Despite 
these changes, the normalized wavelength is quite invariant, similar to the results in 
Tables 1 and 2. The efficiency coefficient 7 also exhibits httle variabihty. It is relatively 
constant with values in the range 0.7 — 0.8 up to 6% of agar concentration, and only for 
8% it slightly drops. The value of 2% of agar concentration corresponds to a standard 
value for which substrate is prepared and data in Tables 1 and 2 were obtained for this 
particular value. 

Relationships between locomotory parameters. 

In order for the worms to move forward (or backward), the wave of neuro-muscular 
activity must travel backward (or forward). One of the main characteristics of this 
wave is its wavelength. In Fig. 4 we plot the normalized wavelength as a function of 
frequency, for a population of wild-type C. elegans, for different mutants, and for related 
Caenorhabditis species. As could be expected from the data in Tables 1-3, these results 
for the forward locomotion show that the normalized wavelength \/L is practically 
frequency independent, and moreover, its value is conserved across all three cases and 
is around 2/3. For backward movement, we observe that X/Lis also about 2/3 (data not 
shown) . This result is consistent with the similarities in the neuroanatomical structure 
of the forward and backward motion neural circuits (White et al, 1986). 

In Fig. 5, we investigate the relationship between the velocity of worm's propulsion 
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and the velocity of the neuro- muscular wave given by XiojI'K. We plot the data for a 
population of wild-type G. elegans, for different mutants, and for related species. We 
find that these two velocities scale linearly with the proportionality coefficient (least- 
square fit slope) around 0.8 across all three conditions. The fact that this coefficient is 
always < 1 implies that the velocity of worm's propagation is less than the velocity of 
the neuromuscular wave, and this is a consequence of the fact that C. elegans worms 
slip during sinusoidal motion. However, as can be seen from Tables 1-3 and Fig. 5, the 
slippage, defined as 1 — 7, is small and to a large extent independent of wavelength, 
frequency, velocity, amplitude, and even the firmness of the substrate. It also does not 
change during development. 

Developmental data provide a broader parameter space for worm's body length and 
other related parameters. In particular, we find that as worm's length increases, the 
amplitude and wavelength increase in a correlated manner. They are related by a linear 
scaling (Fig. 6), which implies that during development their ratio is conserved, i.e., 
Aq/X — 1/5. However, for adult mutants and species, there is some variability in this 
relationship, because the amplitude depends also on frequency and other biophysical 
parameters (see below). 

A particular class of una- 54 mutants studied by us have significantly reduced fre- 
quency of undulations as compared to wild-type worm, and yet they move with a slightly 
greater amplitude of the wave (Table 2; and for an extreme example, see Fig. 1, panel 
d) . These data suggest that there may be some inverse relationship between the ampli- 
tude and frequency of undulations. To test this empirical observation, we determined 
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the statistical significance of this inverse relationship by computing the correlation be- 
tween the amplitude and frequency for the forward motion for the wild-type C. elegans 
larvae and adults and on different substrate, for mutants, and related species (Tables 4 
and 5). We computed correlation coefficients for both instantaneous values of amplitude 
and frequency and their average values. Indeed, these two parameters are negatively 
correlated in all cases and for almost all worms with the same genotype. There are a few 
exceptions, but because the p-values in those cases are large, they should be excluded 
as not statistically significant. In general, correlations for instantaneous values of am- 
plitude and frequency are smaller but more statistically significant than those for their 
average values. The lower negative correlations and lower p-values for instantaneous 
values probably reflect much larger sample size and consequently smaller variability in 
instantaneous locomotory parameters than in averaged parameters. Negative correla- 
tions are also present for worms with different genotypes: among different mutants of 
C. e/e^'ans including wild-type the correlation coefficient is equal to —0.293, and among 
different Caenorhabditis species the correlation is of comparable value —0.329. Similar 
negative correlations are present also for combined forward and backward motion (data 
not shown). 

Theoretical results and their relation to experimental data 

What could be a possible explanation of the constancy of the normalized wavelength? 
It was suggested many years ago by L. Byerly and R. Russell (cited in White et al, 1986) 
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that the wavelength may be related to the spatial extent of motor neuron's dendrites 
that have presumably stretch receptors at their endings (Fig. 2). These researchers 
concluded that since the dendrite length is fixed then so should be the wavelength. It can 
be shown, using a coupled oscillator model (Kopell and Ermentrout, 1988; Cohen et al. 
1992; J. Karbowski, unpublished results), that the Byerly- Russell hypothesis provides 
only a partial explanation of this conservation phenomenon. In fact, the wavelength of 
the emerging wave depends not only on the spatial extent of long-range coupling but 
also on the nature of the interaction function, which in turn, depends on the underlying 
biophysics. The latter factor could be potentially frequency dependent (Gutkin et al, 
2005), but only for frequencies that are much larger, of the order of 50 Hz, than relevant 
frequencies for nematode locomotion, which are ~ 1 Hz. This probably can explain the 
independence of the wavelength on frequency for wild- type C. elegans and among related 
species, which have similar nervous systems. The matter is more subtle with C. elegans 
mutants. It is natural to expect that mutations not directly affecting neural function 
should not influence long-range coupling. On the other hand, mutations investigated by 
us which do affect neural function, affect only local inter-neuronal synaptic transmission. 
These mutations do not affect dendritic structure, and thus, presumably do not alter 
activities of stretch receptors. It thus follows that all our mutations should not influence 
long-range coupling between body segments, consistent with the conservation of the 
normalized wavelength. 

The experimental flndings from Fig. 5 impose two constraints on a theory: linear 
scaling between the velocities of propulsion and neuromuscular wave, and the constancy 
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of the slope. The mechanical part of our model enables us to gain insight about the 
conditions under which the two empirical constraints are met. Computation of a worm's 
propulsion velocity v yields a following relationship between v and the neuromuscular 
velocity, Aa;/27r (see, Appendix): 



.^T^. (13) 



where the propeller efficiency coefficient 7 is given by 



_ [{cn/cl) - 1] (27rAo/A)^ 
^ 2 + {cN/cL){27rAo/X)^ ' 



The result for 7 is essentially the same as that derived for flagellar swimming on the 
basis of the Resistive Force theory (Childress, 1981). The fact that v is linearly related 
to uj follows from neglecting worm's mass and associated inertial forces. The coefficient 
7 depends on the ratio of frictional coefficients in the normal Cat and longitudinal cl 
directions of the main body axis, and on the ratio of the amplitude and wavelength. In 
the case relevant for worm's locomotion, 7 is positive and in the range between and 
1. (7 is negative for cl > and in this case worm moves in the same direction as does 
neuromuscular wave, which contradicts empirical observations and is not considered by 
us.) The value of 7 approaches its maximal value 1, if any of the above ratios becomes 
large. Since the amplitude is weakly frequency dependent and can change by a factor 
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of two among different mutants and species (Table 2), one could expect that the coef- 
ficient 7 may be frequency dependent, and consequently, the relationship between the 
velocities of propulsion and muscular wave may have a non-linear character. However, 
the experimental data show that these two quantities scale linearly with almost the 
same coefficient 7, close to the optimal value 1, across wild- type and many mutant C. 
elegans, and across related species (Fig. 5). Moreover, the ratio 2t:Aq/\ is of the order 
of 1 (Fig. 6, and Tables 1-3). This experimental hnear scaling can be reconciled with 
the formula (13), only if the ratio cn/cl of the frictional coefficients is large, such that 7 
becomes a frequency-independent constant approaching its maximal value 1. Since the 
experimental data give 7 ~ 0.8 for almost all cases, we can estimate from this the ratio 
of the frictional coefficients between worm's body and a substrate as cjq /cl ~ 9.0 — 14.0. 
Surprisingly, this ratio is preserved even on substrates with increased firmness up to a 
large agar concentrations, which is evident from conserved 7 in those cases (Table 3). 
This conservation presumably takes place because Cat and change proportionally on 
these substrates such that their ratio remains constant. 

A few mutants have shghtly lower values of the coefficient 7. These are sqt-1 mu- 
tants, defective in a cuticle collagen, and cat-4 mutants, which arc depleted in dopamine 
and serotonin and have weakened cuticle (Loer et al, 1999) (C. Loer, pers. comm.). 
They have the proportionality coefficient 7, around 0.5 — 0.6, which corresponds to the 
significantly reduced ratio ciq/cL ~ 3.0 — 5.9. This reduction is presumably due to some 
structural changes in the cuticle arrangement that impact its external condition, which 
in turn affect friction. However, even these exceptions still preserve the proportionality 
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relation between both velocities within their respective populations. 

From Eq. (6), we determine the magnitude and spatio-temporal dependence of the 
total bending moment M{x, t). It depends non-linearly on the position along the worm's 
length taking values zero at both endings and it exhibits wave-like property (Fig. 7; 
and Appendix) . This wave corresponds to the neuro-muscular wave traveling from head 
{x — 0) to tail {x — L). Moreover, the magnitude of M is proportional to the difference 
in the neuro-muscular and propulsion velocities, i.e., it is proportional to 1 — 7. For 7 
close to its maximal value 1, which is the case for C. elegans, its mutants, and related 
species, the magnitude of the total bending moment approaches zero (see Appendix). 
Thus, the efficient /optimal locomotion corresponds to the disappearance of the total 
bending moment M and the shear force F. In this limit, the muscle bending moment 
exactly balances body reaction represented by the visco-elastic bending moment 
Me, and worms do not slip. 

Experimental data show negative correlations between the amplitude and frequency 
of undulations (Tables 4 and 5). We seek to understand this inverse relationship in 
theoretical terms. To achieve this, we introduced explicitly simplified muscle dynamics 
into our model of undulatory movement and coupled it with the mechanics. This 
coupling allows us to determine the dependence of the amplitude of undulations on the 
frequency and other parameters (see Appendix): 



^0 = 



A 



KfT 




(14) 
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Thus, the amphtude depends on several macroscopic and microscopic parameters: the 
visco-elastic properties of the hydrostatic skeleton and cuticle characterized by E and 
II, the neuromuscular synaptic transmission T, the rates characterizing contraction- 
relaxation cycle of muscle activity / and g, and the muscle stiffness k. The correction 
O in Eq. (14) is of the order of (1 — 7), which is small for the efficiency coefficient 
7 is close to its upper limit 1. In general, formula (14) implies that the amplitude is 
a decreasing function of the frequency given that all other parameters are constant, 
which is the case for the worms with identical genomes. However, this dependence 
can be weak if characteristic frequencies are smaller than parameters characterizing 
relaxation processes E/iJ, (in cuticle/hydrostatic skeleton) and f + g (in muscles). This 
can explain weak negative correlations between the amplitude and frequency among 
worms representing the same genotype. Our experimental data are consistent with 
the decaying trend predicted by Eq. (14) both among worms representing the same 
species and among worms representing the same mutation (Fig. 8), although detailed 
fits are not possible due to large noise in the behavioral data. This noise arises from the 
changes in the worm's direction of movement. The fact that negative correlations are 
also present among different mutants and among different species representing different 
genomes (Tab. 4) indicates that the parameters k, /, g, T, E, /i, and J in Eq. (14), do 
not vary dramatically in these worms. The possible variability of these parameters is 
not strong enough to reverse the decaying trend of the amplitude with the frequency. 

Prom our theoretical analysis it follows that the amplitude in Eq. (14) is proportional 
to the wavelength A. Thus, if frequency of undulations and other parameters in this 
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equation do not change much, then the amphtude should scale linearly with A. 
Developmental locomotory data provide a good opportunity to verify this scaling, since 
A changes by a factor of 4 during development. Indeed, data from Fig. 6 confirm our 
theoretical result. 

Experimental data from Table 2 indicate that the set of hyperactive mutants with 
affected synaptic signaling {goa-1, egl-30, cat-2, and partly cat-4) move faster and with 
higher frequencies than wild-type worms. These behavior is opposite to a more common 
behavior exhibited by most other neuronal mutations which cause locomotory slowdown. 
Such decreased velocity and frequency is also typical for mutants with structural changes 
in cuticle (BE109 and sqt-1) and muscle {unc-54)- Although our biomechanical model 
cannot explain why frequency changes, it can explain the cross- dependencies between 
locomotory parameters. We suspect that frequency increases in hyperactive neuronal 
mutants because they have elevated overall synaptic transmission that activates neu- 
rons and muscles more vigorously. On the other hand, for sluggish neuronal mutants 
frequency drops presumably because synaptic transmission is less effective. Muscle and 
cuticle structural mutants move with lower frequency probably because they have mod- 
ified neural activity via mechanosensory feedback. Because velocity and amplitude of 
undulations depend on frequency, once the frequency is altered, they change accord- 
ingly. In the case of velocity, there exists its linear relationship with frequency (Fig. 
5), and all mutants conform to that simple rule. The matter is more complex with the 
amplitude dependence on frequency, because several other factors are also involved (see 
Eq. 14). This complexity likely explains why the general trend of negative correlations 
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within a given genotype (Tables 4 and 5, and Fig. 8) does not always translate directly 
to such a dependence across different genotypes (Table 2). For example, neuronal mu- 
tants in comparison to wild- type, despite having higher frequencies, do not have smaller 
amplitude. These mutants presumably have increased overall synaptic transmission T 
at the neuromuscular junction (see Eq. (14)), which counterbalance and even overcomes 
the effect associated with the frequency increase. Structural cuticle and muscle mutants 
have more predictable amplitude in comparison to the wild-type worms. Cuticle mu- 
tants have either comparable or increased amplitude, because their frequency is lower 
and their elasticity coefficient E is slightly reduced, causing Aq to increase in relation 
to wild- type (see Eq. 14). Similarly, muscle unc-54 mutants have in general increased 
amplitude because they have reduced frequency and because of the inverse dependence 
of Aq on the crossbridge dissociation rate g. In unc-54 mutants this rate is presum- 
ably decreased in comparison to wild- type (Moerman and Fire, 1997). Some unc-54 
dominant-negative mutants have altered muscle structure and are paralyzed (Bejsovec 
and Anderson, 1988, 1990). These mutants may have severely reduced muscle stiffness 
coefficient k, which leads to ~ and consequently the propulsion velocity v ~ 
(see Eq. (13)). Since the amplitude of undulations is proportional to the synaptic 
transmisson T at the neuromuscular junction and to the muscle contraction rate /, so 
is the efficiency coefficient 7. However, increase in either of these parameters will not 
necessarily lead to a significant increase in 7 and v, because 7 is already close to a 
saturation caused by a large ratio of cn/cl- Increase in 7 and velocity would be much 
more pronounced in environments in which cn and cl are comparable. 
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General discussion and summary 

We sought to systematize experimental behavioral data by providing a coherent 
model that could help us in understanding the main characteristics of the nematode 
movement and their interdependencies. In particular, our model points out the causes 
of the sinusoidal crawling and hints about its stability 

Our major experimental findings are: (i) the velocity of worm's propulsion scales 
linearly with the velocity of the neuro-muscular wave with almost optimal and highly 
conserved proportionality coefficient, (ii) the value of the normalized wavelength is con- 
served across a population of C. elegans, their mutants, and across related Caenorhab- 
ditis species, and (iii) the amplitude of undulations is inversely correlated with the 
frequency of the wave. These data provided enough information to constrain our model 
and to derive some conclusions about parameters controling the undulatory motion. 

Our experimental data indicate that the frequency of undulations varies between 
different Caenorhabditis species, between different mutants of C. elegans, and even be- 
tween different worms representing the same genotype. This observation suggests that 
the presumed oscillatory activity of motor neurons can be easily modulated, which can 
be explained by noting that nerve ring input and presumably mechanoscnsory receptors 
are involved in modulating ocillatory activity of motor neurons. Since both of these 
contributions can depend on the inner state of an animal, it is not surprising that the 
frequency of undulations can be altered for the worms with the same genotype. For ex- 
ample, we expect that mutations to genes involved in inter-neuronal transmission {goa-1 
and egl-30 encoding proteins that regulate synaptic transmission, cat-2 involved in syn- 
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thesis of dopamine, and cat-4 involved in synthesis of dopamine and serotonin) affect the 
input coming to motor neurons and consequently modulate the frequency of their activ- 
ity. However, our biomechanical model cannot explain precisely why frequency changes 
under different conditions. This would require an approach that explicitly takes neural 
dynamics aspects into account. Nevertheless, our model can explain relationships of 
other locomotory parameters to the frequency. 

Nematode worms do not possess limbs and therefore have to use other strategies for 
efficient locomotion. These worms move primarily because of two main factors: genera- 
tion of the neuro-muscular wave, and the presence of oscillations in local units between 
activities of dorsal and ventral parts of the body. The second factor, characterized by 
the frequency, is the driving force for locomotion, since it enables the wave to acquire 
speed (equal to Xuj/2'k) and consequently to propagate along the body. The traveling 
wave is crucial for the movement, since if it does not propagate, then the worm does 
not move (non-propagating wave has Xuj = and, from Eq. (13), v = 0). In the 
case when there are no oscillations, the neuromuscular wave can still exist (provided 
stretch receptors are active) but only as a standing wave that bends the body, which 
however is not enough to generate worm's movement. Indeed, C. elegans often pause, 
and maintain their body posture, then resume movement. 

The wave can propagate in either direction of the body and this is strictly related 
to the direction of worm's motion. If the wave propagates backward, then the worm's 
muscles exert force on the ground that is directed toward the tail. The mechanics 
principle of action and contraction (the environment reacts with the opposite force on 
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the worm's body) is responsible for the worm's forward motion. If the neuro-muscular 
wave propagates forward along the body, then the directions of the forces are reversed, 
and consequently the worm moves backward. Prom neurophysiological data it is evident 
that the direction of the wave is correlated with the direction of dendritic processes of 
motor neurons in the forward and backward motion neural circuits (Chalfie et al, 1985; 
White et al, 1986). This observation was the motivation for Byerly and Russell (cited 
in White et al, 1986) to propose that dendrites may contain stretch receptors that are 
used in wave generation. 

Oscillations in local units create anti-phase activity between dorsal and ventral parts 
of the body, which in turn, is responsible for a nonzero amplitude of undulations. The 
nonzero value of this amplitude is necessary for transforming part of the ncuro-muscular 
wave velocity into the velocity of propulsion. This can be seen from Eq. (13), which 
relates the two velocities and contains the explicit dependence of the propeller efficiency 
coefficient 7 on the undulatory amplitude. The greater the amplitude of undulations, 
the greater the efficiency of movement (i.e., the propulsion velocity is greater for a given 
velocity of the neuro-muscular wave). However, this increasing trend has its limit, since 
amplitude and frequency are negatively correlated (see Tables 4 and 5). Thus, at higher 
frequencies the amplitude decreases, and the movement efficiency would decrease. The 
fact that the slopes of the regression lines and some data points in Fig. 4 are close to 
the maximal allowed value for 7 suggests that Caenorhabditis worms are evolutionarily 
adapted to maximize this coefficient. However, based on our estimation of parameters, 
they do it not by increasing the amplitude but by increasing the ratio of the frictional 
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coefficients cn /cl- This presumably takes place because of the longitudinal arrangement 
of structural elements in the cuticle that influences its external condition such that it 
enables worms to move much easier longitudinally than normally to the main body 
axis. This choice of the optimizing parameter is more beneficial, since it potentially 
maximizes 7 in different environments regardless of the values of other parameters. 

Although all mutations examined by us change the average frequency of undulations, 
they do not have disastrous effects on movement. (There are of course mutations that 
abolish movement; Brenner, 1974.) This observation suggests that the stability of 
the neuro-muscular traveling wave is robust against perturbations, and our model can 
provide an explanation for it. To destroy the traveling wave, one would have to either 
eliminate oscillations or to alter the long-range coupling between worm's body segments 
which is presumably mediated by long neuronal processes. Our genetic perturbations 
do not change the anatomical structure of the nervous system, instead, some of them 
change local chemical signaling between neurons. This, however, only modulates the 
wave speed by changing its frequency, and does not influence the wave stability. 

We found that the amplitude and frequency of undulations are negatively corre- 
lated. Especially, the correlation between instantaneous values of these parameters is 
important, since it captures sudden changes in worm's movement. All wild-type species 
and most C. elegans mutants show robust negative instantaneous correlations with very 
small p- values (p < 10^"^). However, for some unc-54 mutants p- values are large, which 
suggests that correlations are not statistically significant in those cases. This can be ex- 
plained by the formula in Eq. (14). unc-54 worms move much slower than the wild-type 
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and their frequency is small. Equation (14) implies that for sufficiently small frequen- 
cies, amplitude of undulations becomes practically frequency independent. Thus, any 
sudden change in frequency would have almost no effect on amplitude, which explains 
why these two parameters are poorly correlated in these cases. 

The mechanical part of our model assumes that components of the resistive force 
acting on nematodes are proportional to the corresponding components of the velocity 
of motion. This assumption, called the Resistive Force Theory, is borrowed from fluid 
dynamics models describing swimming (Gray and Hancock, 1955; Lighthill, 1976). In 
our experiments, the worms do not swim. Instead they crawl through a thin layer 
of bacteria and therefore the law of the Resistive Force Theory may be vahd only 
approximately. It would be interesting to investigate a more realistic resistance law 
applicable to our locomotory system, although experimentally it may be hard to achieve 
due to microscopic sizes of nematodes. 

In summary, we presented an integrated approach to studying the undulatory loco- 
motion of nematode worms. Since physiological approaches are difficult in C. elegans, 
our experimental approach consisted of genetic perturbations to parameters controling 
the movement. By comparing different mutants and Caenorhabditis species we were 
able to construct a biomechanical theoretical model that provides insight about differ- 
ent factors involved. Our results reveal that although different mechanical parameters 
characterizing the undulatory locomotion change as we apply perturbations, their in- 
terdependencies are robust and do not fall apart easily. In particular, we have found 
optimahty of the motion reflected in the value of the efficiency coefficient close to its 
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maximal value. This optimality corresponds to almost exact balance between the muscle 
bending moment and the visco-elastic body reaction moment. Conserved relationships 
between locomotory parameters suggest robust cellular and molecular controlling mech- 
anisms that can be fleshed out with the extensive knowledge of details of the C. elegans 
nervous system. 



Appendix 

In this Appendix we provide more details on some derivations. Solution of Eq. (6) 
in the main text yield the following spatio-temporal dependence of the total bending 
moment M{x, t): 



M{x^ t) — ^ — (1 — 7) [sin(a;i — kx) — sinlcut) + kx cosiut) 

-{kx^/L) [cos{ujt - kL) + 2cos{ujt) + {3/kL){sm{Lot - kL) - sin(cjt))] 
+{kx^/L'^) [cos{u;t - kL) + cos{u;t) + {2/kL){s\u{ujt - kL) - siu{ujt))] 



(15) 



where k — 27r/A. M(x, t) is composed of the traveling wave given by sm.{ijjt—kx) and the 
standing wave (the rest of the terms). The latter originates from the recoil correction 
5h{x, t). It is important to note that the total bending moment is proportional to 1 —7, 
and thus for the efficiency coefficient 7 approaching its upper limit 1, the moment M 
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is zero. 

The explicit form of the visco-elastic bending moment Mg{x, t) obtained by inserting 
Eq. (7) into Eq. (8) is: 



Me(x, t) — — [E cos{ujt — kx) — jJLU sin(a;t — kx)\ . (16) 
A 



This moment as well as the total bending moment are both proportional on the ampli- 
tude Aq of undulations. Its space averaged form is: 



Me{t) = ^I^^^ {Esmujt + imjj cosut) , (17) 



where we used experimental fact that A ~ 2L/3. 



The spatially averaged muscle bending moment is determined from solving Eqs. 
(10)-(12) in the main text. It has the form: 



KfT 



^rn{t) = ,,,, , , , [(/ + g) COs{ut + + ^ + 0] ■ (18) 



The dependence of the amplitude of undulations on the frequency cu can be 
determined from the above equations using the facts that M = Mg + M^, and M pa 
in the limit 7 ~ 1, which is justified experimentally. Thus, in this limit we have a 



perfect balance of and Mg, i.e., Ri — Mg. Because both moments and M, 
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oscillate in time, the next step is to take squares of both sides of this equation and then 
to perform temporal averaging over one period of oscillation, i.e.. 



(Mj) = (Me> + 0(l-7)- (19) 

The averaging is done using the facts that {sin^ cot) = {cos^ ujt) = 1/2 and (sina;^ cosa;^) = 
0. As a result of this procedure, we obtain Eq. (14) in the main text. 

The propulsion velocity v is determined from Eq. (9) using Eqs. (7) and (15) for the 
lateral displacement h and the total bending moment M, respectively. After tedious 
algebra we obtain: 

/•L r)M rP'h 1 

/o a^a^^ = -^cr,AlkL{u - vk) [l + O {l/ikLf)] , (20) 
where the main contribution comes from the traveling wave part of M, and 

r ^"^^f ^ -\^l^kL [1 + {l/{kL))] , (21) 

where again the main contribution comes from the traveling wave part of h{x,t). The 
correction to the velocity from the recoil term is small, of the order of l/{kL) 0.1. In 
general, v is approaching the velocity of the neuromuscular wave, Aa;/27r, Hcn/cl ^ oo. 
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Tabic 1. 



Forward locomotory data for wild-type C. elegans during different developmental stages. 



Development 


N 


V [mm/s] 


cu/27r [Hz] 


Aq [mm] 


A [mm] 


X/L 


7 


stage 
















LI 


4 


0.04±0.02 


0.26±0.12 


0.04±0.01 


0.18±0.02 


0.71±0.05 


0.93±0.61 


L2 


4 


0.08±0.01 


0.34±0.06 


0.06±0.00 


0.27±0.02 


0.59±0.03 


0.82±0.18 


L3 


5 


0.08±0.02 


0.26±0.06 


0.08±0.01 


0.36±0.01 


0.63±0.01 


0.84±0.28 


L4 


5 


0.15±0.05 


0.31±0.12 


O.lliO.Ol 


0.54±0.01 


0.66±0.02 


0.89±0.45 


Adult 


5 


0.12±0.03 


0.20±0.07 


0.14±0.03 


0.70±0.03 


0.65±0.02 


0.86±0.37 



N is the number of worms used in every developmental stage. Parameter 7 is defined 

as 7 = f / (Aa;/27r). 
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Tabic 2. 

Forward locomotory data for adult wild-type (WT) C. elegans, its mutants, 
and related species. 



Genotype N v 


mm/s] u!/2tt [Hz] Aq 


mm 


A 


mm 


A/L 7 



a elegans WT 


58 





17±0 


04 





36±0 


08 





09±0 


01 





59±0 


04 





62±0 


02 





79±0.26 


Mutants: 






































cat-2(elll2) 


13 





24±0 


04 





44±0 


06 





11±0 


01 





67±0 


04 





62±0 


02 





82±0.18 


cat-4(ell41) 


10 





12±0 


03 





43±0 


06 





12±0 


02 





57±0 


02 





62±0 


02 





48±0.13 


egl-30(tg26) 


16 





19±0 


03 





57±0 


08 





10±0 


00 





41±0 


02 





54±0 


02 





81±0.17 


goa-l(nll34) 


12 





22±0 


05 





53±0 


10 





09±0 


01 





54±0 


04 





63±0 


02 





76±0.23 


goa-l(syl92) 


11 





24±0 


03 





58±0 


06 





11±0 


01 





53±0 


02 





64±0 


01 





76±0.12 


Ion- 1 (el 85) 


5 





21±0 


03 





40±0 


03 





13±0 


01 





67±0 


03 





62±0 


02 





76±0.12 


Ion- 1 (el 85); 






































lon-2(e678) 


5 





15±0 


03 





33±0 


08 





16±0 


02 





60±0 


05 





56±0 


06 





75±0.24 


BE109 


9 





07±0 


03 





18±0 


06 





11±0 


01 





43±0 


04 





57±0 


04 





91±0.48 


sqt-l(sclOl) 


13 





10±0 


03 





26±0 


07 





12±0 


01 





69±0 


04 





63±0 


03 





55±0.22 


sqt-l(scl03) 


5 





09±0 


03 





27±0 


10 





09±0 


02 





53±0 


04 





62±0 


04 





62±0.31 


unc-54(stl30) 


5 





02±0 


01 





05±0 


00 





11±0 


02 





55±0 


04 





69±0 


03 





73±0.36 


unc-54(stl32) 


5 





08±0 


01 





15±0 


03 





10±0 


01 





59±0 


02 





65±0 


03 





90±0.21 


unc-54(stl34) 


5 





05±0 


01 





08±0 


01 





10±0 


01 





56±0 


03 





65±0 


03 


1 


16±0.28 


unc-54(stl35) 


5 





02±0 


01 





05±0 


01 





16±0 


02 





57±0 


01 





62±0 


02 





60±0.32 


unc-54(s95) 


5 





04±0 


01 





07±0 


01 





12±0 


03 





63±0 


04 





66±0 


02 


1 


00±0.29 


unc-54(s74) 


5 





04±0 


01 





08±0 


02 





13±0 


01 





48±0 


02 





58±0 


02 





97±0.34 



Species: 

C. briggsae 5 

C. japonica 5 

C. remanei 5 

PSIOIO Caen. 4 

CB5161 Caen. 5 



0.15±0.06 
0.08±0.03 
0.25±0.06 
0.19±0.01 
0.14±0.03 



0.30±0.12 
0.15±0.04 
0.41±0.09 
0.39±0.04 
0.29±0.05 



0.09±0.02 
0.15±0.02 
0.15±0.02 
O.lliO.Ol 
0.13±0.02 



0.56±0.02 
0.69±0.04 

0.67±0.04 
0.54±0.04 
0.67±0.06 



0.59±0.02 
0.68±0.03 

0.64±0.03 
0.68±0.03 
0.62±0.02 



0.91±0.51 
0.71±0.32 
0.93±0.30 
0.93±0.12 
0.71±0.20 
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Tabic 3. 

Forward locomotory data for adult wild-type C. elegans on different substrates. 



Agar percentage v [mm/s] uj/2tt [Hz] A, 



mm 



A mm 



A/L 



7 



2 % (N=18) 
4 % (N=18) 
6 % (N=17) 
8 % (N=17) 



0.21±0.04 
0.17±0.04 
0.15±0.04 
0.10±0.02 



0.39±0.08 
0.34±0.06 
0.34±0.06 
0.27±0.05 



0.13±0.02 
0.11±0.02 
0.12±0.03 
0.13±0.02 



0.72±0.03 
0.64±0.05 
0.62±0.07 
0.61±0.05 



0.66±0.02 
0.62±0.04 
0.62±0.03 
0.62±0.04 



0.75±0.21 
0.77±0.23 
0.71±0.24 
0.63±0.17 
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Table 4. 



Correlation between the amplitude Aq and frequency uj of for the forward 
motion of wild-type C. elegans, its mutants and related species. 





Correlation coefficient: 


Correlation coefficient: 




instantaneous Aq and u) 


average Aq and u) 


C. elegans Wild-type 










LI stage 


-0.417 (0.000 


3268) 


-0.812 (0.188 


4) 


L2 stage 


-0.169 (0.000 


2410) 


-0.904 (0.096 


4) 


TO J_ 

LS stage 


-0.290 (0.000 


4278) 


-0.770 (0.128 


5) 


L4 stage 


-0.350 (0.000 


4308) 


-0.966 (0.008 


5) 


Adult 


-0.251 (0.000 


68302) 


-0.505 (0.000 


58) 


Mutants: 










cat-2(elll2) 


-0.286 (0.000 


16034) 


-0.714 (0.006 


13) 


cat-4 (el 141) 


-0.149 (0.000 


7872) 


-0.444 (0.199 


10) 


egl-30(tg26) 


-0.206 (0.000 


10431) 


-0.228 (0.500 


11) 


goa-l(nll34j 


-0.300 (0.000 


12922) 


-0.537 (0.072 


12) 


goa-l(syl92) 


-0.244 (0.000 


11017) 


-0.391 (0.235 


11) 


1 -4 / -4 r \ 

lon-l(el85) 


-0.214 (0.000 


17296) 


-0.799 (0.105 


5) 


lon-l(el6oj; lon-2(eo loj 


-0.104 (0.000 


10670) 


0.266 (0.665; 


5) 


1)77 1 r\r\ 


-0.168 (0.000 


10000) 


-0.256 (0.505 


9) 


sqt-l(sclul J 


-0.189 (0.000 


14038) 


-0.442 (0.131 


13) 


sqt-l(SClUoj 


-0.382 (0.000 


5551) 


-0.963 (0.008 


5) 


r 1 / J. -4 (Tin \ 

unc-54(stl3u) 


-0.002 (0.918 


3785) 


-0.935 (0.020 


5) 


unc-54(stl32) 


-0.237 (0.000 


4968) 


-0.027 (0.966 


5) 


unc-o4(Stlj4 J 


-0.024 (0.083 


5457) 


-0.629 (0.256 


5) 


unc-5Afstl35) 


0.0068 (0.675; 3824) 


-0.184 (0.767 


5) 


unc-54(s95) 


-0.026 (0.058; 5506) 


-0.955 (0.011 


5) 


unc-54(s74) 


-0.060 (0.001; 2965) 


-0.915 (0.029 


5) 


mutants cross-correlation 






-0.293 (0.247 


17) 


Species: 










C. briggsae 


-0.432 (0.000; 7061) 


-0.974 (0.005; 5) 


C. japonica 


-0.255 (0.000; 4659) 


-0.790 (0.112; 5) 


C. remanei 


-0.380 (0.000; 4955) 


-0.335 (0.582; 5) 


PSIOIO Caenorhahditis 


-0.497 (0.000; 5049) 


-0.831 (0.169; 4) 


CB5161 Caenorhahditis 


-0.266 (0.000; 5172) 


0.236 (0.703; 5) 


species cross-correlation 






-0.329 (0.262; 6) 
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For each correlation coefficient we show in the bracket corresponding 
p-value and the number of data pair points, respectively. 
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Table 5. 



Correlation between the amplitude Aq and frequency uj of for the forward 
motion of adult wild-type C. elegans on different substrates. 



Agar percentage 


Correlation coefficient: 


Correlation coefficient: 




instantaneous Aq and u 


average Ao and uj 


2% 


-0.325 (0.000; 17378) 


-0.246 (0.325; 18) 


4% 


-0.213 (0.000; 19421) 


-0.099 (0.697; 18) 


6% 


-0.183 (0.000; 17828) 


-0.270 (0.294; 17) 


8% 


-0.155 (0.000; 17950) 


-0.383 (0.129; 17) 



For each correlation coefficient we show in the bracket corresponding 
p-value and the number of data pair points, respectively. 
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Figure Captions 
Fig. 1 

Photographs of typical shapes of worms performing undulatory locomotion recorded in 
our experiments. An extreme example is unc-54(stl35) worm, which moves very slowly, 
and with low frequency but has a large amplitude of the wave (panel d). The white 
midline with dots on it was used as a reference frame for quantifying the movement. 

Fig. 2 

The system controling locomotory output in Caenorhabditis. Nerve ring neurons acti- 
vate motor neurons, which in turn activate body-wall muscles. Muscle activity causes 
the worm to move and the movement is modulated by hydrostatic skeleton and cuticle, 
and by mechanosensory feedback. The latter component provides stability for the un- 
dulations. The lower panel shows basic components of the neural structure for forward 
locomotion: the nerve ring, ventral nerve cord (dashed line), and excitatory B motor 
neurons. Neuronal processes in B motor neurons (only 3 shown) are elongated and it 
is hypothesized that their endings contain stretch receptors. These long dendrites pro- 
vide long-range coupling between remote segments of the body, and their directionality 
correlates with the direction of motion. In the circuit controling backward motion den- 
drites in A neurons are elongated in the opposite direction. Ventral side of the worm 
corresponds to the side containing cell bodies. 

Fig. 3 

Diagrams showing worm's body and forces acting on a slice of its body during undu- 
latory locomotion. All the forces and moments acting on the slice must balance each 
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other. 



Fig. 4 

Conservation of the normahzed wavelength \/L across a population of wild- type C. 
elegans, its mutants, and across related Caenorhabditis species for forward movement. 
Diamonds represent data points and solid lines are the least-square fits to the for- 
mula X/L = Oo + a-iuj. (A) For a population of wild type C. elegans the fit yields: 
Qq = 0.629, Oi = —0.035 {N = 58). Each data point corresponds to average values 
oi X/L and cu for one worm. (B) For different mutants of C. elegans including its 
wild-type the fit yields: oq = 0.639, ai = —0.068 (N = 17). The mutants included: 
BElOl sqt-l(sclOl) and BE103 sqt-l(scl03), BE109, BC347 unc-54(s74), RW130 unc- 
54(stl30), RW132 unc-54(stm), RW134 unc-54(stl34), RW135 unc-54(stl35), and 
RW5008 unc-54(s95), goa-l(nll34) and goa-l(syl92), egl-30(tg26), cat-2(elll2), cat- 
4(ell41), CB185 lon-l(el85), and PS697 lon-l(el85);lon-2(e678). The data points are 
averages over a population for each mutation. (C) For related Caenorhabditis species 
the fit yields: Oq = 0.659, Oi = —0.066 [N = The species included: N2 C. elegans, 
AF16 C. briggsae, SB339 C. japonica, SB146 C. remanei, CB5161 Caenorhabditis sp., 
and PSIOIO Caenorhabditis sp. The data points are averages over a population for each 
species. Note that the normalized wavelength is almost frequency independent (oi is 
close to zero) and very similar in all three figures. 

Fig. 5 

Linear scaling of the velocity of forward propulsion (velocity of the center of mass) with 
the velocity of muscle contraction wave and conservation of the coefficient 7. Diamonds 
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represent data points and solid lines are the least-square fits to Eq. (13) in the main 
text. The slope of the regression line corresponds to the average coefficient 7. The slope 
of the dashed hnes indicate the maximal allowed slope, i.e., 1. (A) For a population of 
wild type C. elegans: 7 = 0.787 (A^ = 58, i?^ = 0.900). Each data point corresponds 
to average values for one worm. (B) For a group of mutants of C. elegans, the same as 
in Fig. 4B, including its wild-type: 7 = 0.791 {N = 14, R'^ = 0.988), except for 3 data 
point (circles) representing BElOl sqt-l(sclOl), BE103 sqt-l(scl03), and cat-4(ell41) 
mutants that have significantly reduced the ratio cn/cl- (C) For related Caenorhabditis 
species, the same as in Fig. 4C, the slope is 7 = 0.861 {N — Q,B? — 0.882). The data 
points in (B) and (C) are averages within a population for each mutation (B) and each 
species (C). Note a very similar and almost optimal value of the coefficient 7 for all 
three cases. 

Fig. 6 

Linear scaling of the amplitude of undulations with the wavelength during different 
developmental stages. The least-square fit to the data points yields regression line 
= 0.194A + 0.006 with i?^ = o.98. 

Fig. 7 

Dependence of the total bending moment M on the position along worm's body. This 
dependence has a non- linear character with a travelling wave of activity present. Pa- 
rameters used: Aq = O.l mm, L = 1.0 mm, A = 0.66 mm, uj = 1.8 Hz, 7 = 0.8, Cat = 50 
g/ (mm s). 
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Fig. 8 

Dependence of the amplitude of undulations on the frequency for different species (A), 
(B), and C. elegans mutants (C), (D), (E). Diamonds are data points and solid lines are 
the least-square fits to the formula Aq — a/ [{1 + 6^a;^)(l + c^ou'^)]^^^ , which is equivalent 
to that represented by Eq. (14). One data point corresponds to one worm. The 
parameters of the fits are: (A) for AF16 C. briggsae yields a — 0.121, b — 2.903, 
c = 0.099 (A^ = 5, i?2 = 0.934); (B) for SB339 C. japonica a = 0.201, b = 6.207, 
c = 0.001 {N = 5, = 0.560); (C) for BE103 sqt-l(scl03) mutants a = 0.140, 
b = 3.977, c = 1.466 (A^ = 10, R'^ = 0.830); (D) for BC347 unc-54(s74) mutants 
a = 0.155, b = 9.478, c = 0.040 (A^ = 5, i?^ = 0.871); (E) for cat-2(elll2) mutants 
a = 0.209, b = 3.673, c = 0.004 {N = 13, i?^ = 0.503). The fit for BC347 unc-54(s74) 
allows us to estimate an effective crossbridge dissociation rate, which turns out to be 
about 2.5 times smaller than in the other two C. elegans mutants: sqt-l(scl03) and 
cat-2(elll2), which do not alter muscle dynamics. 
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